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Abstract. The problem of building supersymmetry in the quantum mechanics of two 
Coulombian centers of force is analyzed. It is shown that there are essentially two ways of 
proceeding. The spectral problems of the SUSY (scalar) Hamiltonians are quite similar and 
become tantamount to solving entangled families of Razavy and Whittaker-Hill equations 
in the first approach. When the two centers have the same strength, the Whittaker-Hill 
equations reduce to Mathieu equations. In the second approach, the spectral problems are 
much more difficult to solve but one can still find the zero-energy ground states. 

Key words: supersymmetry; integrability; quantum mechanics; two Coulombian centers 
2000 Mathematics Subject Classification: 81Q60; 81V55 



1 Introduction 

Supersymmetry is a bold idea which arose in the seventies in string and field theory. It was 
immediately realized that mechanisms of spontaneous supersymmetry breaking should be inves- 
tigated searching for explanations of the apparent lack of supersymmetry in nature. In a series 
of papers [1, 2, 3] Witten proposed to analyze this phenomenon in the simplest possible setting: 
supersymmetric quantum mechanics. A new area of research in quantum mechanics was born 
with far-reaching consequences both in mathematics and physics. 

Of course, there were antecedents in ordinary quantum mechanics [nihil novum sub sole), and 
indeed even before. The track can be followed back to some work by Clifford on the Laplacian 
operator, see [4], quoted in [5]. Recast in modern SUSY language, the Clifford supercharge is: 

Q = ( n m+ "I" + ^ ), 2 2 = j 2 = k 2 = — 1, ij = -ji = k ciclyc, 



n 
where i, j, k are the imaginary unit quaternions and 

are the components of the gradient modified by the components of the electromagnetic vector 
potential. The SUSY Hamiltonian is: 

otn + not - ( ~ A + iB ^) + 3 B ^) + k Mx) 
^V-hW y Q - A —iB\{x) — jB 2 (x) — kB 3 (x) J ' 



*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
http:/ /www. emis.de/journals/SIGMA/symmetry2007. html 
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the Laplacian plus Pauli terms. Needless to say, an identical construction relates the Dirac 
operator in electromagnetic and/or gravitational fields backgrounds with the Klein-Gordon 
operator. The factorization method of identifying the spectra of Schrodinger operators, see [6] 
for a review, is another antecedent of supersymmetric quantum mechanics that can also be 
traced back to the 19th century through the Darboux theorem. 

In its modern version, supersymmetric quantum mechanics prompted the study of many 
one-dimensional systems from a physical point of view. A good deal of this work can be found 
in [7, 8, 9, 10]. A frequent starting point in this framework is the following problem: given 
a non-SUSY one-dimensional quantum Hamiltonian, is it possible to build a supersymmetric 
extension? The answer to this question is positive when one finds a solution to the Riccati 
equation 

ldWdW hd^W_ _ 
2 dx dx 2 dx 2 

that identifies the - a priori unknown - "superpotential" W(x) from the - given - potential V{x). 
Several examples of this strategy have been worked out in [11]. 

The formalism of physical supersymmetric systems with more than one degree of freedom 
was first developed by Andrianov, Ioffe and coworkers in a series of papers [12, 13], published in 
the eighties. The same authors, almost simultaneously, considered higher than one-dimensional 
SUSY quantum mechanics from the point of view of the factorization of TV-dimensional quan- 
tum systems [14, 15]. Factorability, even though essential in TV-dimensional SUSY quantum 
mechanics, is not so effective as compared with the one-dimensional situation. Some degree of 
separability is also necessary to achieve analytical results. For this reason we started a program 
of research in the two-dimensional supersymmetric classical mechanics of Liouville systems [16]; 
i.e., those separable in elliptic, polar, parabolic, or Cartesian coordinates, see the papers [17] 
and [18]. We followed this path in the quantum domain for Type I Liouville models in [19], 
whereas Ioffe et al. also studied the interplay between supersymmetry and integrability in quan- 
tum and classical settings in other type of models in [20, 21]. In these papers, a new structure 
was introduced: second-order supercharges provided intertwined scalar Hamiltonians even in the 
two-dimensional case, see [23] for a review. This higher-order SUSY algebra allows for new forms 
of non- conventional separability in two dimensions. There are two possibilities: (1) a similarity 
transformation performs separation of variables in the supercharges and some eigenfunctions 
(partial solvability) can be found, see [24, 25]. (2) One of the two intertwined Hamiltonian 
allows for exact separability: the spectrum of the other is consequently known [26, 27]. 

The second difficulty with the jump in dimensions is the identification of the superpotential. 
Instead of the Riccati equation one must solve the PDE: 

\vW(x) • VW{x) + ^V 2 W{x) = V(x). 

Z Z 

In our case, we look for solutions of this PDE when V{x) is the potential energy of the two 
Coulombian centers. We do not know how to solve it in general, but two different strategies 
should help us. First, following the work in [38] and [39] on the supersymmetric Coulomb 
problem, we shall choose the superpotential as the solution of the Poisson equation: 

^V 2 W{x) =V(x). 
z 

The superpotential will be the solution of another Riccati-like PDE where a classical piece must 
be added to the potential of the two centers. Second, as in [41, 34] the selection of superpotential 
requires the solution of the Hamilton-Jacobi equation: 

^VW(x)- VW(x) = V(x). 
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Again the superpotential must solve a third Riccati-like PDE, where now a quantum piece must 
be added to the potential of the two centers. 

The organization of the paper is as follows: We start by briefly recalling the non-SUSY 
classical, Newtonian [40], and quantal, Coulombian [42, 33], two-center problem. We shall 
constrain the particle to move in one plane containing the two centers. The third coordinate is 
cyclic and it would be easy to extend our results to three dimensions. In Section 2 the formalism 
of two-dimensional SUSY quantum mechanics is developed, and the superpotential of the first 
Type is identified for two Coulombian centers. Bosonic zero-energy ground states are also found. 
Sections 3 and 4 are devoted to formulating the SUSY system in elliptic coordinates where the 
problem is separable in order to find fermionic zero-energy ground states. It is also shown that 
the spectral problem is tantamount to families of two ODE's of Razavy [43], and Whit taker-Hill 
type [44], see also [28]. Since these systems are quasi-exactly solvable, several eigenvalues are 
found following the work in [32]. In Section 5 two centers of the same strength are studied and 
some eigenfunctions are also found. Section 6 is fully devoted to the analysis of the Manton- 
Heumann approach applied to the two Coulombian centers. Finally, a summary is offered in 
Section 7. 

1.1 The classical problem of two Newtonian/Coulombian centers 




Figure 1. Location of the two centers and distances to the particle from the centers. 

The classical action for a system of a light particle moving in a plane around two heavy bodies 
which are sources of static Newtonian/Coulombian forces is: 



rn 



dx\ dx\ dx2 dx2 \ ot\ 0^2 



dt dt 



+ 



dt dt J 



T2 



The centers are located at the points (x\ — — d, #2 = 0), ( x i = X2 = 0), their strengths are 
a\ = a > o>2 = 5a > 0, 6 G (0, 1], and 



r i = Y ( x i - °0 2 + x h r 2 = y ( x i + °0 2 + x 2 

are the distances from the particle to the centers. In the following formulas we show the dimen- 
sions of the coupling constants and parameters and define non-dimensional variables: 



[ai] = [a 2 ] = [a] = ML 3 T~ 2 , 
X\ — > dx\, X2 — > dx2, t 



[d] 



[5] 



d 3 m 



a 



ri — > dr\ = d\j (x\ - l) 2 + x 2 , r 2 — > dr 2 = dJ {x\ + l) 2 + x\. 



In the rest of the paper we shall use non-dimensional variables. From the non-dimensional action 

dx2 dx2 \ 1 5 
ri r 2 



S — V mda S — V mda 



( 1 / dx\ dx\ 
\2 \ ~dt~dt 



+ 



dt dt 
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the linear momenta and Hamiltonian are defined: 
dL dx\ dL dx2 



Pi 



dx\ dt 



P2 



3x2 dt 



2 T\ T2 



This system is completely integrable because there exists a "second invariant" in involution with 
the Hamiltonian: 



h (mda)I 2 , h = ~{l 2 ~ pi) + %i 



S_ 

T2 



1 



l 2 = (xip 2 - x 2 pi) 



1.2 The quantum problem of two Coulombian centers of force 



If \fmad is of the order of the Planck constant ft, the system is of quantum nature. Canonical 
quantization in terms of the non-dimensional ft constant, 

d 



Pi Pi 



-ih 



dxi 



[xi,pj] = ihS, 



ft = 



X{ > X{ 

h 



%ii 



\JmdoL 

converts the dynamical variables into operators. The quantum Hamiltonian, H — , and the 
quantum symmetry operator, I2 — (mda)l2, are mutually commuting operators: 

[Hj 2 \ = HI 2 -i 2 H = Q, 



* h ( d 2 d 2 

H = ~2\dxi + 



dx\ 



1 6 

+ — + — . 
n r 2 



l2 = ~Y 



d 2 



dx\ X2 dx\ 



- 2X\X2 



d 2 



d 



d 



a a ~ Xl x X2 ^~ ) + Xl 

OX\OX2 OX\ OX2 



r 2 n 



2 Two-dimensional Af = 2 SUSY quantum mechanics 

We now describe how to build a non-specific system in two-dimensional J\f = 2 SUSY quantum 
mechanics. Besides commuting - non-commuting - operators there are anti-commuting - non- 
anti-commuting operators to be referred respectively as "bosonic" and "fermionic" by analogy 
with QFT. The Fermi operators are represented on Euclidean spinors in R 4 by the Hermit ian 
4x4 gamma matrices: 
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{ 7 J ,7 fc } = W k = {7 2+J ,7 2+fc }, W,l 2+h } = 0, j,k = 1,2. 
The building blocks of the SUSY system are the two (A/* = 2) quantum Hermitian super- 
charges: Q\ = Q u Ql = Q 2 , 



2+j 



(y)t y 5 



(7 



2+j\t 



4 - v^t 



dW 

^ ~ ^ 2 



2 ( 
.7=1 ^ 



8x, 



J=l J= 
It is convenient to define the non- Hermitian supercharges Q± = Q± ± zQ2, 
/ \ 









dXj 



dW 
dxi 



fc a dW 

8X2 0X2 



V 



_k_d_ 1 <9VK t, d 



9X1 u / 
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dW 



o 

















_t_a aw 

8x2 dx2 



9xi 







because their anti-commutator determines the Hamiltonian H$ of the supersymmetric system: 

{Q+,Q.} = 2hH s , [Q+,H s ] = [Q-,Hs]=0. 
The explicit form of the quantum SUSY Hamiltonian is enlightened by the "Fermi" number 

2 

F — Yl ^h-V 7 - operator: 



/ 



fcg> fcg 



/i 



(1) 

21 



V o 









(i; 

12 

(1) 
22 





o \ 



U 2 ) j 



F 



i=i 



/ \ 

10 

10 

V 2 y 



It has a block diagonal structure acting on the sub-spaces of the Hilbert space of Fermi numbers 
0, 1, and 2. In the sub-spaces of Fermi numbers even Hg acts as ordinary differential Schrodinger 
operators. The scalar Hamiltonians are: 

2/>(/=°) = -n 2 v 2 + vwvw + nv 2 w = -h 2 v 2 + 2V {0 \ 
2h {f=2) = -h 2 v 2 + vwvw - hv 2 w = -h 2 v 2 + 2V {2) . 



In the sub-space of Fermi number 1, however, Hg is a matrix of differential operators, the 2x2 
matrix Hamiltonian: 



2h«=V 



-n 2 v 2 + vwvw - nu 2 w 



-2h 



d 2 W 



-2h 



V = ei + e 2 , 

OX\ OX2 



a 2 w 

dx\dx2 



v 2 = 



dx\dx2 

-Wv 2 + vwvw + m 2 w 
d 2 d 2 , d 2 



+ 



dx\dx\ 8x28x2' 



□ 2 = 



d 2 



dx\dx\ 8x28x2 



This is exactly the structure unveiled in [12, 13]. All the interactions expressed in H$ come from 
the as yet unspecified function W(xi,X2), which is thus called the superpotential. 



2.1 The superpotential I for the two-center problem 

To build a supersymmetric system containing the interactions due to two Coulombian centers of 
force, we must start by identifying the superpotential. One possible choice, inspired by [38], is 
having the two-center potential energy in in the term proportional to h. We must therefore 
solve the Poisson equation to find the superpotential I: 

2 n T2 h h 

Note that the anticommutator between the supercharges induces a H factor in front of the 
Laplacian. This fact, in turns, forces the singularity of the superpotential (henceforth, also of 
the potential) at the classical limit h = 0. The same singularity arises in the bound state spectra 
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Figure 2. Cross section (x 2 = 0) and 3D graphics of the quantum potential for S = 1/2. Cases: 
Upper row: (a) ^ = 0.2, (b) h = 0.4. Lower row: (a) h = 1 and (b) ^ = 10. Increasing ft the centers 
become more and more attractive. 



of atoms, e.g., in the energy levels of the hydrogen atom. The potential energies in the scalar 
sectors are accordingly: 



v^ 2)) = 



2 

¥ 



1 + S 2 + 6 



-) 

r 2 n nr 2 l 



n T2 



l s 

T | - + - 
n r 2 



We stress that the superpotential I is a solution of the Riccati-like PDE's: 

viyviy ± w 2 w = 2vy (2); 

and the scalar Hamiltonians for two SUSY Coulombian centers read: 



W(0)\ 



£ 2 . 



-V 2 + 



fi 2 



1 + S 2 + 5 



-) 

r 2 ri rir 2 / 













T2J 



2.2 Bosonic zero modes I 

The bosonic zero modes 

0±*o 0) (^2) = 0, Q T y ( o\x 1 ,x 2 ) = 0, 
if normalizable, are the bosonic ground states of the system: 

j exp[ (- 2 ^ 2) ] \ / 



V 









V$ t \x 1 ,x 2 ) 



I 



V exp[ 







(2ri+2fr 2 ) i 
k 2 J 



The norm of the true bosonic ground state \I/q^ is finite and given in terms of Bessel functions: 

jw( X1 , X2 ) = 2 r 

Jo 

4 , „A „ / 4 



/ dxidx 2 en w{ - Xl ' X2) =2 I dx x \ dx 2 e^ 

-oo </-oo </0 JO 

w 



2 (-2n-25r 2 ) 



^- ar li(rTfl*b (1 "' ) r I ^ (1+ ' ) 

+ 4(i^(!< i - i >M! (1+i> 

In Fig. 3 plots of the zero energy bosonic ground state probability density of finding the 
particle in some area of the plane are shown for several values of ft. The drawings reveal the 
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physical meaning of the h = singularity: for h = 0.2 we see the particle probability density 

peaked around the center on the right with a very small probability. Exactly at h — 0, e~V 
is finite only for T\ = and zero otherwise, giving probability of finding the particle exactly in 

the center. e _ ~^~, however, is zero W2 meaning that the probability of this state is zero at the 
classical limit; in classical mechanics there are no isolated discrete states. 




Figure 3. Graphics of the norm as a function of h and the ground state probability density, 
|^"q°^ (xi , 0^2) | 2 , for S = 1/2 and the values of h = 0.2, 0.4, 1 and 10. Note the extreme smallness 
for h = 0.2. The norms for these four cases are: 7V(0.2) = 3.5806 • 1(T 47 , N(0A) = 2.0576 • 10 -13 , 
N(l) = 0.00942 and iV(10) = 1743.94. 



3 Two-dimensional Af = 2 SUSY quantum mechanics 
in elliptic coordinates 

The search for more eigenfunctions of the SUSY Hamiltonian requires the use of the separability 
in elliptic coordinates of the problem at hand. This, in turn, needs the translation of our two- 
dimensional Af — 2 SUSY system to elliptic coordinates. A general reference (in Russian) where 
SUSY quantum mechanics is formulated in curvilinear coordinates is [36], see also [37] to find 
a more geometric version of SUSY QM on Riemannian manifolds. 
The change from Cartesian to elliptic coordinates, 

x\ — uv G (—00, +00), X2 — =b\/ (u 2 — 1)(1 — v 2 ) E (—00, +00), 
u=^(ri+ r 2 ) e (1, +00), v = ^(r 2 - n) e (-1, 1), 

induces a map from the plane to the infinite elliptic strip: M? = (— 00, +00) x (—00, +00) =^> 
E 2 = (—1, 1) x (1, +00). This map also induces a non-Euclidean (but flat) metric in E 2 : 

g{u,v) = 

with Christoffel symbols: 




—v 



= -u(l -v 2 ) = v(u 2 - 1) = 

uu (u 2 - v 2 )(u 2 - 1)' vv (u 2 - v 2 )(l - v 2 )' uv vu u 2 -v 2 ' 
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uu 



v(l-v 2 ) 



VV 



-u(u 2 



uv 



(u 2 - v 2 ){u 2 - 1) ' vv [u 2 - V 2 ){1 - V 2 ) ' 

Using the zweig-bein chosen in this form, 

2 2 

g uu (u, v) = e ]( u > v)e](u, v), g vv (u, v) = ^ e)(u, v)e v j (u, v), 

3=1 3=1 



vu 



U 2 — V 2 



ei(u,v) 



u - 1 

u 2 — V 2 



e%(u,v) 



1-v 2 
u 2 — V 2 



we now define "elliptic" spinors, "elliptic" Fermi operators, and "elliptic" supercharges: 

i>l(u,v) = et{u,v)ip l ±, if>l(u,v) = e v 2 (u,v)^ 2 ± , 

( 0\ 



C+ = -iVh 



C- = -iVh 



















V 



v ofa+- J^) \ 



hu 



(0 ^(V+ + 





\0 



hu 



v u 2 —v 2 







where W(u, v) = F(u) + 

To obtain the supercharges in Cartesian coordinates from the supercharges in elliptic coor- 
dinates, besides expressing u and v in terms of x\ and x%, one needs to act by conjugation with 
the idempotent, Hermitian matrix: 



S 



( 1 \ 

— ve\{u, v) —ue^iujv) 

— ue^u, v) ve\(u,v) 

\ o o o -l / 



SC+S = Q+, SCS = Q- 



Equation (1) in elliptic coordinates, 
■ 2 1 (d 2 F u dF^ 



u 



+ 



u 2 — v 2 \ du 2 u 2 — 1 du 



+ 



d 2 G 



u 2 — v 2 \ dv 2 



v dG 
— v 2 dv 



(1 + S)u (8-l)v 



u 2 — V 2 



+ 



u 2 — V 2 



(2) 



is separable: 



2 d 2 F ^^dF + 2(1 + S)u 
du 2 du h 



(1 - v 2 ) d2 G + v d ^ + 2 ( 5 ~ l ) v 
dv 2 dv ft 



with separation constant n. The general solution of equation (2) depends on two integration 
constants (besides an unimportant additive constant): 



W(u, v; «, C u C 2 ) = -2^r^u + £ In (u + V^ 2 - l) + ^ fin (u + ^ 2 - 1 

a a 2a V 

+ 2 



, C 2 . K 2 

- — v H — - arcsm v = ( arcsm v) . 

h h 2h 
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We find thus a three-parametric family of supersymmetric models for which the potentials in 
the scalar sectors are: 



~ ( ( 2 } ) 

Vj (2) (xi,x 2 ;k, C\,C 2 ) 



i + s 2 + s(^ + r ^- 

T2 T\ 



-) 



+ 



Ci + n In 



n + r 2 + \J (ri + r 2 ) 2 4 



+ a 



ft arcsm 



r 2 - T\ 



2(1 + 5)V(n + r 2 ) 2 -4 Ci + ft In 



n + r 2 + v (ri + r 2 ) 2 - 4 



+ 2(1 - 5)^4 - (r 2 - ri) 2 C 2 - ftarcsin 



r 2 - ri 



1 5 

T I - + - 

ri r 2 

We shall restrict ourselves in the sequel (as before) to the simplest choice ft = C\ = C 2 = such 
that we shall work with the "elliptic" superpotential I 



W(u,v) 



2(1 + S)u 2(1 -6)v 



h 



h 



W(x 1 ,x 2 ) 



2ri 25r 2 



because this election is significative and contains enough complexity. 

Considering families of superpotentials related to the same physical system in our 2D frame- 
work differs from a similar analysis on the ID SUSY oscillator, see e.g. [29], in two aspects: 
(a) Because we solve the Poisson equation, not the Riccati equation, the family of superpo- 
tentials induces different families of potentials in both and h^ 2 \ (b) and are not 
iso-spectral because are not directly intertwined. The spectrum of is the union of the spectra 
of &°) and ti 2 \ 



3.1 Fermionic zero modes I 

The fermionic zero modes 

( \ 



/ o \ 

V / 



( \ 
V o J 



A, 



y/u 2 — v 2 



e n 




+ 



A, 



are fermionic ground states if normalizable. Because the norm is: 



/l roo 
dv J 



du 



V(« 2 -i)(i-« 2 ) 



> F(u) — G(v) 
' h + 



\Ju 2 — V 2 



^( U 2_1)(1_ W 2) 



/ o \ 



F{u)-G{v) 

e n 

V o J 



n F(u)-G(v) 

e z n 



it is finite if either A\ — or A 2 = 0. With our choice of sign in F(u) the fermionic ground 
state is: 



\/u 2 — V 2 






2(l+<5)u+2(l-<5)v 
P 



\ 
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and the norm is also given in terms of Bessel functions: 



N(h) 



roc rl 



dv 



e n 



Vu 2 - 1 v / T 3 " 



2^o(|(l + 5))/ (|(l-5) 



It is also possible to give the fermionic ground state in R 2 using the § matrix: 

/ \ 

_2(Sr 1 ±r 2 ) 



4 _ ri _ r2 _j_ 2 



\ 



(ri+r 2 ) 
4 v / rir2 

; (r2-ri) 

T2 V rir2 T2 ri 





fi2 




Figure 4. Graphics of N(H) for: (5 = 1/2, « = 0. ^(zi, x 2 )\ 2 for (5 = 1/2, H = 0.2, 1, 4 and 10. 
Norms: 7V(0.2) = 1.3518 • 10~ 45 , N(l) = 0.0178, iV(4) = 7.3881, iV(10) = 18.4297. 



4 The bosonic spectral problem I 

The spectral problem for the scalar Hamiltonians is also separable in elliptic coordinates. Plug- 
ging in the separation ansatz 



( ( 2 } ) ( ( 2 } ) ( ( 2 } ) ( ( 2 } ) 



in the above spectral equation we find: 

-h 2 (u 2 - 1)^ - ^ 2 «^ + (4^|^V - 1) T 2(1 + 5)« - 2Eu 2 

((0)\ 

" ^ + h*v± + (4^2(1 - v 2 ) T 2(1 - 5), + 2Ev 2 

((0)\ 

= -4 (2) >), 

where / is the eigenvalue of the symmetry operator / = -{/i'( 2 )' + 7 2 (2) }. 



((0)\ 
H2)// \ 



((0)\ 

d (2) >) 
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Research on the solution of these two ODE's by power series expansions will be published 
elsewhere. Here, we shall describe how another change of variables transmutes the first ODE 
into Razavy equation [43, 28]: 

d 2 ri±(x) . ^ _ _ x2 / \ /\ 1 

— ^ V ((± cosh2x — M±) T]±{x) A± ?7±(£), x -arccosh^ 

ax z 2 

with parameters: 



M± 



2 2(1 + 5f 



2(1 + S) 2 - K 2 E± 



Simili modo, another change of variables leads from the second ODE to the Whittaker-Hill or 
Razavy trigonometric [44, 28], equation 

+ (/?±cos2y - N±f Z±(y) = fi± $±(y), y = l - arccost; £ [0, f ] 
with parameters: 



fc-^a-*-^, ni- w-o- 



h\ h 2K J ± 2(1 - 5) 2 -h 2 E ± 



h 2 V ^ 2 

Therefore, the spectral problem in the scalar sectors is tantamount to the solving of two 
entangled sets - one per each pair (E,I) - of Razavy and Whittaker-Hill equations. If M± = 
nf + 1, nf E N + , the Razavy equation is QES; i.e., there are known n + 1 finite eigenfunctions 
with an eigenvalue, see [32]: 

This means that for those values of E± one expects bound eigenstates of the SUSY Hamil- 
tonian, although the ^-dependence cannot be identified. If N± = + 1, E N + , there exist 
finite eigenfunctions of the Whittaker-Hill equation, with eigenvalues: 

Again, one expects these values of E± to be eigenvalues of the SUSY Hamiltonian, although 
their eigenfunctions are expected to be non-normalizable (except the ri^ — case) because the 
^-dependent part of the eigenfunction is non finite and u is a non-compact variable. In any 
gives the bosonic zero mode. 



4.1 The fermionic spectrum I 



The eigenfunctions of the matrix Hamiltonian, 
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l 21 



(i) _ *(i) 

12 

(1) 



(xi - l) 2 - x\ ^ [x\ + l) 2 - x\ 
x 2 {x\ - 1) x 2 {x\ + \) 



-2 



+ 8- 



i 



22 
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except the fermionic ground states, follows easily from the SUSY algebra. The fermionic eigen- 
functions in elliptic coordinates, 



V(u,v) = C + ^> + (u,v) 
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in Cartesian coordinates {x\,x<i) = S^^ ± (u,v) read: 
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5 Two centers of the same strength 

If the two centers have the same strength, 5 = 1, the spectral problem in the scalar sectors 
becomes tantamount to the two ODE's: 



-*V - - ^ + (1(- 2 - 1) T 4n - 2£n 2 



d 

"du 
d 



h 2 



((0)\ /(0)\ 



.H\l-v 2 )^ + h 2 v^ + 2Ev 2 
dv z av 



((0)\ /(0)\ 



Identical changes of variables as those performed in the < 6 < 1 cases now lead to the Razavy 
and Mathieu equations [35], with parameters: 

~ rf2 ^2 X) + (C± cosh 2s -M±) 2 r/± 0*0 = \± v ±(x), 
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c± = ± ; 
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d 2 Z±(y) 
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X± = Mi + w [I ± + 
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+ (a±cos4y + a±) £±(y) = 0, a± = 



4E± 



We now select the three lowest energy levels from finite solutions of the Razavy equation: 



Eq = 0, Ei 
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The corresponding eigenfunctions of the Razavy Hamiltonians for ri^ = 0, 1, 2 are: 
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The energy degeneracy is labeled by the eigenvalues of the symmetry operator 



/r = | (Ar _ (n + 1)2) _16 ) 



m = 1, 2, . . . , n + 1, 
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In this rotationally non-invariant system the symmetry operator / replaces the orbital angular 
momentum in providing a basis of common eigenfunctions with H in each degenerate in energy 
sub-space of the Hilbert space in such a way that a quantum number reminiscent of the orbital 
angular momentum arises. 

Next we next consider even/odd in v solutions of the Mathieu equations 

&» = | (C [aT, ?n, arccos(*;)] + C [af, 9n , arccos(-*;)] ) 

+ [a£™ 9n, arccos(i;)] + S [a n ^ , 9n , arccos(-*;)] ) , 

= y (C [ar,«n, arccos(i;)] - C [a£* , ?n, arccos(-^)] ) 
+ y (5 [a^ m , g nj arccos(»] - S [aj m , g n , arccos(-i;)] ) . 

The reason is that the invariance of the problem under r\ ^ r2 = v ^ — v - the exchange 
between the two centers - forces even or odd eigenfunctions in v. The parameters of the Mathieu 
equations determined by the spectral problem are: 



nm 
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E n + IT 1 



En 



4 h 2 ' Hn 8 2W 

To fit in with the parameters of the Razavy Hamiltonians set by nf = 0, 1, 2 we must choose: 



qo = 0, 4 1 = 0, 
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Therefore, 

.,(0) 



nm even/odd 



n > 0, 



l,2,...,m+l, 



is a set of bound states of non-zero energy of the scalar Hamiltonian of two SUSY Coulombian 
centers of the same strength. The paired fermionic eigenstates are obtained through the action 
of the appropriate supercharge. Since every non-zero-energy state come in bosonic-fermionic 
pairs, the criterion for spontaneous supersymmetry breaking is the existence of a Fermi-Bose 
pair ground state of positive energy connected one with each other by one of the supercharges. 
Our SUSY Hamiltonian has both bosonic and fermionic zero modes as single ground states; 
consequently, supersymmetry is not spontaneously broken in this system. 

In the next figures we show several graphics of bosonic SUSY eigenfunctions for the choice: 
C2 = = 0, ci = d\ = 1. 




Figure 5. Graphics with (ft = 1): \^ even (x 1 ,x 2 )\ 2 ; |^+ e ven(^i^2)| 2 and |^+ dd( x i> x 2)| 2 ; 

l^+even(^1^2)| 2 and |^+odd(^l ? • 



5.1 Bosonic and fermionic ground states for two centers of the same strength 

For comparison, we also plot the bosonic and fermionic ground states for several values of ft. 




Figure 6. Graphics of the probability density |^q (x±, ^2)| 2 for 5 = 1, and the values of ft = 0.2, 1, 4 
and 10. 
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Figure 7. Graphics of ^^(xi, x 2 )\ 2 for 5 = 1, and H = 0.2, 1, 4 and 10. The norms are: N(Q.2) = 
7.7012 • 1(T 88 , N(l) = 0.0009, N(4) = 5.8083, iV(10) = 16.6347. 



6 The superpotential II for the two-center problem 



There is another possibility to find the potential energy of two Coulombian centers in a 2D 
Af = 2 SUSY Hamiltonian. The superpotential must satisfy the Hamilton-Jacobi equation 
rather than the Poisson equation (1): 



1 5 _ 1 (dW\ 2 1 (dW\ 2 



t'i + r 2 2 V dxi J + 2 I dx 2 ) 



(3) 



Note that in this case the two centers must be repulsive to guarantee a real W . This point 
of view, which follows the path shown in [41] and [34] for the Coulomb problem, amounts to 
the quantization of a classical supersymmetric system, only semi-positive definite for repulsive 
potentials. 

Again using elliptic coordinates, the Hamilton-Jacobi equation separates 



V - 1) 



dF 
du 



+ 2(1 + %, K = (l- v 2 )[^-\ + 2(S-l)v 



dG 



V 



by plugging in (3) the ansatz: 

W(u, v; k) = F a (u; k) + G^v] ft), a, b = 0, 1. 
The quadratures 



F a (u-K) = (-1) 



■jf 



du. 



G h {v-K) = {-1) 



/V 
•1 



" x /2(l-«5)u + . 



VT 



show that the separation constant k is constrained in order to find real F a (u) and Gb(v): 
2(1 — 5) < k < 2(1 + 5). Note that there are two different possibilities: a = b and a ^ b. 
A different global sign only exchanges the (0) with the (2) and the (1)1 with the (1)2 sectors. 
The superpotential II is thus given in terms of incomplete and complete elliptic integrals of the 
first and second type, see [30, 31]: 
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The scalar Hamiltonians, both in elliptic and Cartesian coordinates, read: 
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6.1 Type Ha and Type lib two-center SUSY quantum mechanics 

We now present the graphics of the scalar potential for a — 6, a system that we shall call 
Type Ha Af — 2 SUSY two Coulombian centers. By the same token, the system arising from 
the superpotential I will be called Type I Af = 2 SUSY two Coulombian centers. 




Figure 8. 3D graphics of the quantum potential for a = b = 1 (or for a = b = 0). We choose 
S = 1/2, n = 3. Cases: (a) h = 0.2, h = 2. Observe that when h is increased this provides a reduction 
of the strength of the repulsive centers and the left center becomes attractive whereas the right center is 
still repulsive, (b) h = 4 and ft = 10, both centers are attractive. With increasing h the centers become 
more and more attractive. 

If a ^ 6, we shall call the system Type lib Af = 2 SUSY two Coulombian centers. The scalar 
potential is drawn in the graphics below for several values of H. 

The differences between the Type Ha and Type lib potentials, increasing with h, are shown 
in the next figure. 

Below we shall discuss the coincidences and differences of the three distinct points of view. 



6.2 The spectral problem 

Both for the Type Ila and Type lib systems the spectral problem in the scalar sectors is separable 
in elliptic coordinates: 

h^^> (u, t,) = E^> (u, v) , 4 (2)j v) = (tOGT (v), 
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Figure 9. 3D graphics of the quantum potential V"(°) for a = 1, b = (or for a = 0, b = 1). 

We choose S = 1/2, n = 3. Cases: (a) ft = 0.2, h = 2. Observe that once again increasing ft provides 
a reduction of the strength of the repulsive centers and the left center becomes attractive whereas the 
right center is still repulsive, (b) ft = 4 and ft = 10, both centers are attractive. 





Figure 10. 3D graphics of the quantum potential for a = 1 = b and a = 1, b = in the cases 
h = 1, h = 10. In this range we note the differences between Type Ha and lib quantum potentials. 
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The separated ODE's are, however, much more difficult (non-linear) than in the Type I system 
and there is no hope of finding explicit eigenvalues and eigenfunctions. 



6.3 Bosonic ground states 

We shall therefore concentrate on searching for the ground states. First, the bosonic zero modes: 

The separation ansatz 

r(o)\ /(o)\ /(o>\ 

4 (2 >V^rWo (2) >) 
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makes these equations equivalent to: 
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The solutions, i.e. the bosonic zero modes, are: 
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It is not possible to calculate analytically the norm in these cases, but we offer a numerical 
integration of the F = bosonic ground states: 
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exp 
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in the next figures. We observe that there is one normalizable bosonic ground state of zero 
energy of Type Ha and one of Type lib once the value of a is set to be one. Remarkably, the 
Type lib zero mode disappears (the norm becomes infinity) at the classical limit. 



10000 
8000 
6000 
4000 



Figure 11. Numerical plot of iV(ft; 3, 1, 1) as Figure 12. Numerical plot of iV(ft; 3, 1, 0) as 
function of ft for S = 1/2. function of ft for 5 = 1/2. 

Some 3D plots of the Type Ha and Type Ha bosonic zero modes for several values of ft are 
shown in the next figures: 




Figure 13. 3D graphics of the ground state probability density \^°\xi, ^2)| 2 , for S = 1/2, k = 3, and 
a = b = 1 (or \^( 2 \x 1 ,x 2 )\ 2 for a = b = 0). Cases: ft 0.2, ft = 2, ft = 4 and ft 10. The norms are 
7V(0.2) = 0.004914, N(2) = 2.83912, JV(4) = 22.8914 and JV(10) = 511.092. 
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Figure 14. 3D graphics of the ground state probability density \^°\xi, #2)| 2 > for S — 1/2, k = 3, and 
a = 1, b = (or |^ (2) 0i, x 2 )\ 2 for a = 0, b = 1). Cases: ft = 0.2, ft = 2, ft = 4 and ft = 10. And the 
norms are iV(0.2) = 9.61622 • 10 20 , N(2) = 473.903, N(4) = 287.687 and N(10) = 1399.71. 



6.4 Fermionic ground states 

Second, the fermionic ground states of zero energy: 

C + ^\u,v) = 0, C-*$-\u,v) = 0. 

Unlike the bosonic case where the logic is or instead of and, note that both equations must be 
satisfied by the fermionic zero modes. The separation ansatz 

makes these equations tantamount to 
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The fermionic zero modes have the form of the linear combination: 

/ \ / 
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Because the norm is 



dudv 



roc ri 



F a (u;K,)-G b (v,K) F a (u;K)-G b (v,K) 

h + Aie z h 



only A\ — or A 2 — are normalizable, depending on the choice of a and b. 

Again, these integrals cannot be computed analytically but the outcome of numerical calcu- 
lations is shown in the next figures for several values of ft. 

In this case we see that the fermionic zero mode of Type Ha does not have a classical limit 
whereas the Type lib fermionic ground state behaves smoothly near ft = 0. 

3D plots of fermionic zero modes, of both Type Ha and lib, are shown in the last figures. 
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Figure 15. Numerical plots of N(H; 3, 1, 1) and N(H; 3, 1, 0) for A 1 = 1,A 2 =0 and 5 = 1/2 (N(H; 3, 0, 0) 
and N(H; 3^0, 1) for A x = 0, A 2 = 1). 




/ 



Figure 16. Graphics of \V§ )W (x u x 2 )\ 2 for A x = 1, A 2 = 0, 5 = 1/2, « = 3, and the sign combination: 
a = b = 1. Cases: h = 0.2,2,4 and 10. These graphics also represent I^q 1 ^ 2 ^!, x 2 )\ with A\ = 0, 
A 2 = 1, 5 = 1/2, « = 3, and a = b = 0. The norms are 7V(0.2) = 1.03792 • 10 22 , N(2) = 251.908, 



7V(4) = 45.7495 and iV(10) = 25.207. 





Figure 17. Graphics of the function x 2 )\ 2 for A\ = 1, A 2 = 0, 5 = 1/2, « = 3, and the sign 

combination: a = 1, b = 0. Cases: ^ = 0.2,2,4 and 10. These graphics also represent I^q 1 ^ 2 ^!, x 2 )\ 2 
for Ai = 0, A 2 = 1 and a = 0, b = 1. The norms are 7V(0.2) = 0.05046, N(2) = 1.46657, 7V(4) = 3.62192 
and iV(10) = 9.20207. 



7 Summary 

In this paper we have built and studied two types of supersymmetric quantum mechanical 
systems starting from two Coulombian centers of force. Our theoretical analysis could be of 
interest in molecular physics seeking supersymmetric spectra closely related to the spectra of 
homonuclear or heteronuclear diatomic molecular ions, e.g., the hydrogen molecular ion or the 
same system with a proton replaced by a deuteron. In for instance, the first case, the value 
of the non-dimensional quantization parameter is checked to be ft = 0.7: 

h = 1.05 • 10" 34 kg • m 2 • s"\ Vmda = 1.493 • 10" 34 kg • m 2 • s"\ 

using the international system of units (SI). 

The first type is defined from a superpotential that solves the Poisson equation with the 
potential of the two centers as the source. In this case, the spectral problem is shown to be 
equivalent to entangled families of Razavy and Whittaker-Hill equations. Using the property of 
the quasi-exact solvability of these systems, many eigenvalues of the SUSY system corresponding 
to bound states have been identified when the strengths of the two centers are different. If the 
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strengths are equal things become easier and some bound states are also found. In summary, for 
our simplest choice of Type I super potential the main features of the spectrum are the following: 

• There is an infinite set of discrete energy eigenvalues in the F = Bose sub-space of the 
Hilbert space: 

The ionization energy is: E^ = 2 n ^(l + 5) 2 , the threshold of the continuous spectrum. 

• There is a sub-space of dimension n + 1 of degenerate eigenfunctions with energy E n . The 
eigenvalues of the symmetry operator / 

ti 2 

Inm "J ((Km ~ (n + l) 2 ) , 771 = 1, 2, . . . , 72 + 1, 

where Km are the n + 1 roots of the polynomial that solves the Razavy equation, label 
a basis of eigefunctions in each energy sub-space. 

• In the case 5=1 the eigenfunctions can be explicitly found. Moreover, the system enjoys 
a discrete symmetry under center exchange: v — v (r\ <-> 7*2). The energy eigenfunctions 
come in even/odd pairs of functions of v with respect to this reflection. This fact suggests 
that the purely bosonic two fixed centers Hamiltonian 

H = V 2 + — + — + C 

2m T\ T2 

2 

enjoys a hidden supersymmetry (if the constant C is greater than ( & m ^r) of the kind 
recently unveiled in [45] for classy one-dimensional systems. This hidden supersymmetry 
is spontaneously broken because the even/odd ground states have positive energy: E± — 
C-8^ > 0. 

• There are eigenfunctions in the F = 1 Fermi sub-space for the same eigenvalues E n , n > 0. 
Analytically, the Fermi eigenfunctions are obtained from the Bose eigenfunctions through 
the action of Q+. 

The second type starts from superpotentials solving the Hamilton-Jacobi equation. There 
are two non-equivalent sign combinations giving two sub-classes. Both Type Ila and Type lib 
superpotentials are defined in terms of incomplete and complete elliptic integrals of the first and 
second kind. The superpotential in this approach is no more than the Hamilton's characteristic 
function for zero energy and flipped potential. The separability of the HJ equation in elliptic 
coordinates means that we can find a "complete" solution of this equation. The spectral problem 
is, however, hopeless for this Type. 

All the zero-energy ground states, bosonic and fermionic, Type I and Type Ila/IIb, different 
strengths and equal strengths, have been obtained. The cross section (x2 — 0) of the probability 
density of some of the ground states are shown, for the sake of comparison, in these final Tables. 
It is remarkable that, despite of being analytically very different, Type I and Type II zero modes 
show similar patterns. 
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Table 1. ft = 1, 5 = 1/2. 
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Table 2. ft = 2, 5 = 1/2. 
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Table 3. h = 4, 5 = 1/2. 
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